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Abstract 

We revisit phenomenology of the minimal gauge-mediated model. This model 
is motivated from the SUSY CP and flavor problems. A specific feature of this 
model is that tan/3 is naturally large, since the B term in the Higgs potential 
is zero at the messenger scale. This leads to significant SUSY contributions to 
various low-energy observables. We evaluate the anomalous magnetic moment of 
the muon and the branching ratio of B —* X s ~f taking account of recent theoretical 
and experimental developments. We find that the current experimental data prefer 
a low messenger scale (~ 100 TeV) and gluino mass around 1 TeV. We also calculate 
the branching ratios of B — ► X s l + l~, B s — > u + fi~, and B~ — ► t~V, and show that 
these observables are strongly correlated with each other in this model. 



1 Introduction 



Low-energy supersymmetry (SUSY) is a very attractive model of physics beyond the 
standard model (SM). In the minimal supersymmetric standard model (MSSM), however, 
general SUSY breaking masses of squarks and sleptons induce too large FCNC and/or CP 
violation effects in low-energy observables. These SUSY FCNC and CP problems should 
be solved in realistic SUSY breaking models. 

Gauge-mediated SUSY breaking [1, 2, 3, 4] is one of the promising mechanisms to 
describe the SUSY breaking sector in the MSSM. The SUSY breaking is transmitted to 
the MSSM sector through the gauge interaction, which induces the flavor-blind soft SUSY 
breaking masses of squarks and sleptons. In the minimal gauge-mediated model (MGM) 
[5, 6], the trilinear scalar couplings (A terms) and Higgs bilinear coupling (B term) are 
zero at tree level, and they are induced from radiative corrections of the gaugino masses. 
In this case, dangerous SUSY CP phases are absent, and the SUSY CP problem is solved. 

In Refs. [7, 8], phenomenology of the MGM was studied. One of the specific features in 
the MGM is that tan (3 is naturally large, and significant SUSY contributions are expected 
in various low-energy observables. They considered anomalous magnetic moment of the 
muon, the B — > X s 7 decay, and the B — » X s l + l~ decay. It was shown that the deviations 
of B{B — > X s 7) and anomalous magnetic moment of the muon from the SM predictions 
are strongly correlated in the MGM. 

Recently, theoretical calculations in both B(B — > X s ^) and anomalous magnetic mo- 
ment of the muon have been improved. Also the experimental data was updated. Accord- 
ing to Refs. [9, 10], a theoretical estimation of B(B — > X^) at next-to-next-to-leading or- 
der (NNLO) is lower than the experimental world average at a 1.4 a level. In Refs. [11, 12] 
the SM prediction of anomalous magnetic moment of the muon has been updated by using 
the recent e + e~ — > it + it~ data. The new SM prediction is larger than the experimental 
data at a 3.4 a level. While we can not still conclude that the deviations come from the 
new physics, there is a room for additional new physics contributions. 

In this paper, we revisit phenomenology of the MGM taking account of the recent 
developments of B[B — > X s ^) and anomalous magnetic moment of the muon. Also, we 
evaluate B(B -> X s l + l~), B(B S -> and B(B~ -> r~V). The SUSY contributions 
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to all the observables are enhanced when tan/5 is large. We show that the deviations 
from the SM predictions are strongly correlated in the MGM and the recent result of 
B(B — > X s 7) prefers a low messenger scale and gluino mass around 1 TeV. This model 
will be tested by various low-energy experiments in addition to the SUSY direct search 
at LHC. 

2 Minimal Gauge-Mediated Model 

First, let us briefly review the MGM. In the gauge- mediated SUSY breaking, the SUSY 
breaking is mediated to the MSSM sector through the SM gauge interactions. In the 
MGM, N pairs of SU(2)x, doublet and of SU(3)c triplet chiral superfields are introduced 
as messenger fields. The SU(2)^ doublets and of SU(3)c triplets are assumed to be in 
SU(5) 5 and 5*-dimensional multiplets. The messenger fields couple to a singlet chiral 
field by the following superpotential, 

W = XiSQ&i, (1) 

where S is a singlet, $ 2 (3) and $ 2 (3) are SU(2) L -doublet (SU(3)c-triplet) superfields. The 
scalar component of S develops a vacuum expectation value ({S)), which gives super- 
symmetric mass terms to the messenger fields, Mm, = \i{S). Also the non- vanishing 
F-component of S ((F s )) induces the SUSY breaking in the messenger fields. 

The gaugino masses in the MSSM are generated by one-loop diagrams of the messenger 
fields and given by 

M i = N-±kg(x i ) = M i g(x i ), (2) 

where A = (Fs)/{S) and xi = A/Mm,- The loop function g(xi) is given in Ref [13], 
but we take Xi ^ 1 in this paper so that g(xi) ~ I. Hereafter, the parameters at the 
messenger scale M M are denoted by bar. The sfermion mass terms are generated by 
two-loop diagrams and given by 

< ~ 1 (2C3M3 2 + 2C 2 M 2 2 + Vm?) , (3) 

when Xi ^ 1. Here, Cj and Y are the quadratic Casimir and the hypercharge of sfermions, 
respectively. In the gauge-mediated SUSY breaking, the squarks become rather heavy 
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compared with the sleptons and Higgs bosons since the soft SUSY breaking terms are 
proportional to the gauge coupling constants. In the MGM the trilinear scalar couplings 
(A) and the bilinear Higgs coupling (B) are assumed to vanish at the messenger scale, 

A = B = 0. (4) 

Those A and B terms are induced by the gaugino-loop diagrams at the weak scale, 
and the relative phases between the gaugino masses and the A and B terms vanish, 
arg(MjA*)=arg(Mj_B*)=0. Therefore, the dangerous SUSY CP phases are absent, and 
the SUSY CP problem is solved naturally in the MGM. 

With the boundary conditions at the messenger scale, we can evaluate the SUSY 
breaking parameters at the weak scale by solving the renormalization group equations 
(RGEs). For the ji parameter we do not specify its origin, and it is determined by 
imposing the correct electroweak symmetry breaking at the weak scale, as 

f i 2 = \ 2 t A 2 -m 2 H -±m 2 z (l + 5 H ), (5) 

where Xt is the top-quark Yukawa coupling, ran is the common Higgs boson mass, and 
A 2 and Sh represent radiative corrections defined in Ref. [6]. 

The B term is generated radiatively at the weak scale. It has two types of contribu- 
tions: one is a Higgsino-gaugino contribution (Bq) and the other is an effective A-term 
contribution (Ba)- They are destructive to each other, and the magnitude of the B pa- 
rameter is even smaller than one naively expected. We follow the result in Ref. [6] for 
the numerical calculation of the B parameter. Once the B parameter is fixed, tan (3 is 
determined by the minimization condition of the Higgs potential, 

where is the CP-odd Higgs boson mass. 

Notice that tan/3 is not a free parameter in the MGM, but, is predicted from other 
parameters. Since the B parameter is generated radiatively, tan (3 becomes naturally large. 
Also, sign(_B) is opposite to sign(M 2 ) in most parameter regions. Then sign (M 2 /i tan (3) is 
determined to be positive from Eq. (6). Moreover, in the MGM, all the soft SUSY breaking 
parameters are determined by only three input parameters, A, Mm, and N . These are 
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specific features of the MGM and very important when we consider phenomenology of 
the MGM. 



3 Low-Energy Observables 

3.1 Anomalous Magnetic Moment of the Muon 

The current experimental measurement of the anomalous magnetic moment of the muon, 
a n — (9 ~ 2)^/2, was reported by Muon (g-2) Collaboration [14], 

a^ xp = 11659208.0(6.3) x 10~ 10 . (7) 

Recently, the SM prediction of a M has been updated based on the new data of e + e~ — > 
7r + 7r~. The most recent calculation including new evaluation of the hadronic light-by- light 
scattering contribution is [11, 12] 

af 1 = 11659178.5(6.1) x 10~ 10 . (8) 

The discrepancy between the experimental data and the SM prediction is 

5a^ = a* xp - af l = (29.5 ± 8.8) x 10~ 10 . (9) 

The most recent result indicates a 3. 4a deviation, and it suggests new physics contribu- 
tions. 

In the MGM, the contribution to a M is approximately given by 

S<B ^^_^_^ (MM , (10) 

Here we assume that all the SUSY mass parameters to be equal to M S usy- The chargino- 
sneutrino loop diagram dominates over other diagrams. The SUSY contribution is en- 
hanced when tan (3 is large and its sign is determined by sign(M 2 /i). In the MGM, the 
sign of M 2 /xtan/3 is positive, and Sa^ GM > 0. This agrees with the sign of the current 
discrepancy of a M in Eq. (9). 
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3.2 B -> X s 7 Decay 

The experimental world average £>(£? — > X s7 ) by the Heavy Flavor Working Group [15] 
is 

B(B -> X s7 ) £7>1 . 6GeV = (3.55 ± 0.26) x 10" 4 , (11) 

with a photon energy cut E 1 > 1.6 GeV. The SM prediction of B[B — > X s7 ) has been 
estimated by including NNLO contributions in Ref. [9]. Recently, additional perturbative 
corrections related with the photon energy cut has been calculated [10]. Combining those 
results, the most recent theoretical prediction is given by [10] 

B(B -> X s7 )£ 7>1 . 6GcV = (2.98 ± 0.26) x 10~ 4 . (12) 

The new SM prediction becomes lower than the experimental data, and the deviation 
is 1.4 a. While the deviation is not significant, the recent result opens new room for 
additional new physics contributions. 

The B — > X s7 decay is described by the following effective Hamiltonian [9], 

AC 8 

tfir 7 = -^m;e^wo.w. (is) 

V2 i=1 

where \x\, ~ to&. The SUSY contributions are encoded to the Wilson coefficients, Cj. The 
most relevant operator for the B — > X s7 decay is 

7 = -^m b S L a^b R F^. (14) 

The calculation of B(B — > X s7 ) at NNLO is quite complicated. Here we refer Refs. [9, 10] 
to the formula. 

In the MGM, the dominant SUSY contributions to C 7 come from the charged Higgs 
boson, the chargino, and the gluino loop diagrams. It is known that the charged Higgs 
amplitude is always constructive with the SM contribution, while the chargino and the 
gluino loop amplitudes are either constructive or destructive with the SM one. In the 
mass insertion approximation, the pure-Higgsino (C^ ) , the Higgsino-wino (Cj W ), and 
the gluino contributions (Cf) are written as [6] 

Cf ~ -r 6 A^tan/3^/ ch (/i 2 /m|), (15) 

1 m Q 
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- { ^r b UnP-J^[UM^ m l)-U^/rnl)], (16) 
- ^ p^ { 4^r b tan^f gl (Mi/rnl), (17) 

where is an average squark mass and the loop functions / c h and f g \ are defined as 

„ , , 13 — 7a; 3 + 2x(l — x) , . . 

^ oWF- 3(,-l/ lQS ^ (18) 



l + lOx + x 2 3x(l + x) , . 

/*(*) = ~ 2(x _i)4 + 7^Tjr lo g^ ( 19 ) 



The non-holomorphic correction to the bottom-quark Yukawa coupling is not negligible 
in large tan/3 [16]. The correction is represented as r 6 = 1/(1 + e b ), and e b in it is 

e fe = |^M 3 /itan/3/ (m\ L ,m\ R , M 3 2 ) . (20) 

The loop function / in e b is defined as 

xy\og(x/y) +yz\og(y/z) + zx\og(z/x) 

J{X,y,Z) = 7 T7 w 7 • (21) 

(x -y)(y - z)(z - x) 

The flavor-violating parameter (<5ll)23 is induced by the RGE effect and is written in the 
leading logarithmic approximation as 

With the boundary condition A = 0, A t is induce by the RGE effect and proportional 
to the gaugino masses. As discussed in the previous section, sign (M 2 /i tan (3) is predicted 
to be positive in the MGM. In this case sign (A t fi tan (3) is also positive and the sign of the 
pure-Higgsino amplitude is opposite to that of the charged Higgs contribution. Then the 
dominant SUSY contributions can cancel with each other. The Higgsino-wino and gluino 
contributions are also constructive with the pure-Higgsino one. 

In order to calculate B(B -> X s ^) at NNLO in the MGM, SUSY contributions should 
be matched into Cj beyond LO at the weak scale. But this is beyond the scope of this 
paper. For the numerical calculation, we use the LO SUSY corrections to Cj at the weak 
scale while we have included the SM contributions at NNLO. 
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3.3 B -> X a l + r Decay 

B(B — > X s / + Z~) in a low dilepton invariant mass region, 1 < m\ < 6 GeV 2 , was reported 
by Babar, (1.8 ± 0.9) x KT 6 [17] and Belle, (1.5 ± 0.6) x 10~ 6 [18]. The current world 
average is given by 

B(B -> X s /+/-) 1<m , <6GcV2 = (1.6 ± 0.51) x 10- 6 . (23) 

The branching ratio in a higher dilepton invariant mass region is also measured. In a higher 
run region the branching ratio suffers from long-distance contributions of J/ip and ■0' 
resonances, and a large theoretical uncertainty is expected. However, in < m\ < 6 GeV 2 
region, the short- distance contribution dominates, and the branching ratio is sensitive to 
the new physics contribution. Therefore in this paper, we consider B(B — > X s l + l~) in the 
low dilepton region. The SM prediction of B(B — > X s /i + /i~) at NNLO of QCD including 
a QED correction is given by [19] 

B(B - X s /iV)i< m 2<6Gcv2 = (1-59 ± 0.11) x 10~ 6 . (24) 

The SM prediction is consistent with the current experimental data. 
The B — > X s l + l~ decay is described by the effective Hamiltonian 

AC 

h%* +1 - = h£*> - ^v tb v* [CMOM + cmom] , (25) 

where new operators 9 and O w are defined as 

9 = (sZ^b L )(h»l), (26) 
do = (sIlMih^l). (27) 

In the MGM, the Wilson coefficients Cg and Cio are dominated by the SM contributions, 
and the SUSY particles can contribute to B(B — > X s l + l~) mainly through Cj. As a result, 
the SUSY contributions to B{B — > X s ^) and B(B — > X s /i + /i~) are strongly correlated. 
The branching ratio in the low dilepton mass region is written by [19] 

B(B -> X s /i + /i~)i<m2<6GcV 2 = (2.19 - 0.543^10 + 0.0281i? 7 + 0.0153i? 7j Ri 

+ 0.0686^7^8 - 0.870R 7 R 9 - 0.0128i? 8 + 0.00195i? 8j Rio 

- 0.0993^8^9 + 2.84^9 - 0. 107R 9 R 10 + 11.0i? 2 

+ 0.281.R 2 . + 0.00377^ + 1.53R 2 g ) x 10~ 7 . (28) 
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Here R^s are ratios of the total and the SM contributions, 



(29) 



where Cf s are the effective Wilson coefficients at /^vf — m t an d they are the LO contri- 
butions [19]. 

£?(£> — > X s 7) depends mainly on |CV| 2 , and it is insensitive to the sign of C-j. However, 
B(B — > X s l + l~) depends on the sign of CV since there is an interference term between R 7 
and Rg. When the sign of C 7 is opposite to the SM prediction, B[B — > X s l + l~) becomes 
too large compared with the current experimental data, and such a possibility is strongly 
disfavored [20]. 

3.4 B s — > Decay 

The current upper bound on B(B S — > fi + fi~~) is given by the CDF collaboration [21], 



at 95% C.L. In the SM B(B S — > ) is suppressed by the muon mass and quite small. 
The SM prediction is estimated as B(B S -> /x + /Osm = (3.35 ± 0.32) x 10~ 9 [22]. The 
current experimental upper bound is about two orders of magnitude larger than the SM 
prediction, and a clue to new physics may be hidden there. 
The branching ratio of B s — > is written by [23] 



In Ref. [24], it was pointed out that the amplitude of neutral Higgs boson exchange 
diagram is proportional to tan 3 (3 in a large tan (3 region due to non-holomorphic correction 
to the Yukawa coupling. The dominant SUSY contribution is written by 



B(B S -> < 1 x 10 



,-7 



(30) 




C s ^ -C P ~ 



m b m^m 2 t 



tan 3 (3 At/if (m| 




(32) 
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and C' s = {m s /mb)Cs and C' P = —(m s /mb)Cp. The dominant SUSY contribution to 
B(B S — > is proportional to tan 6 /? and can be enhanced by orders of magnitude 

when tan/3^ 10. Since r b < 1 in the MGM, the tan/5 enhancement becomes mild. For 
the numerical calculation we use the complete expression in Refs. [25, 26]. 

3.5 B — > r~V Decay 

The Belle Collaboration reported evidence of B~ — > t~ 17, and the measured branching 
ratio is [27] 

B(B~ -> T - V) = (l^oitistat.^oiUsyst.)) x 10" 4 . (33) 
In the SM the branching ratio is given by 

B{B -^- v ^<^^y M . (34) 

Using \V ub \ = (4.31 ±0.30) x 10~ 3 [28] and f B = 0.216±0.022 GeV [29], the SM prediction 
is B(B~ — > t~v) = (1.6 ± 0.4) x 10 -4 . The current experimental data is consistent with 
the SM prediction. 

In the MGM, the charged Higgs boson exchange contributes to the B~ — > t~V decay. 
The branching ratio is written by 

B(B - ^ T - v) = « (, _ TB _ Sl yj? (, _ (35) 

where the non-holomorphic correction r& is included [30]. The charged Higgs boson am- 
plitude interferes destructively with SM contribution and become significant when tan (3 
is large. The current experimental data gives the constraint on the charged Higgs contri- 
bution. However, a large tan (5 region is still allowed since the cancellation occurs between 
the SM and the charged Higgs contributions. 

4 Numerical Analysis 

In the MGM, there are three input parameters, A, Mm, and N. The numerical results do 
not depend strongly on the difference between M Ms and M Mi , so we assume that M Ma = 
1.3Mm 2 (= M m ) in the following analysis. The gaugino masses are almost determined 
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by A so we will show the numerical result as a function of the gluino mass M 3 . We 
consider the messenger scale as M M /A = 2, 10, 10 2 , 10 3 , and 10 4 . For the quark masses, 
we use m t = 172.5 GeV and m b = 4.7 GeV [28]. With the inputs at the messenger 
scale, we can calculate the SUSY breaking parameters at the weak scale by solving the 
RGE of the MSSM. As discussed in Section 2, tan (3 and /x are determined from the input 
parameters by imposing the correct electroweak symmetry breaking. Once all the SUSY 
breaking parameters at the weak scale are determined, the SUSY contributions to the 
low-energy observables can be evaluated. For convenience, tan (3 and \x are shown in 
Fig. 1 as functions of M 3 for various Mm's. 

First, we discuss the light Higgs boson boson mass, which is shown in Fig. 2 as a 
function of M 3 for various Mm's. We have calculated the Higgs mass using the FeynHiggs 
code [31]. The Higgs boson mass depends mainly on M 3 (i.e. A). From the LEP bound on 
the neutral Higgs boson mass (m h > 114.4 GeV), we can obtain the following lowerbound 
on the M 3 , 

M 3 > 800 ~ 900 (800 ~ 1100) GeV, for N = 1 (2). (36) 

We find that this gives strong bounds on signatures of the MGM. 

It is pointed out in Ref. [6] that the stability of the charge-conserving vacuum gives a 
lowerbound on M 3 in the MGM. When tan (3 is quite large, the lightest stau becomes light. 
In addition, the trilinear coupling of the left- and right-handed staus with the Higgs field, 
which is enhanced by tan (3, may destabilize the charge- conserving vacuum. However, it 
is found from Fig. 1 in Ref. [6] that the lowerbound on M 3 is about 700 ~ 800 GeV even 
for M M ~ 2A when N — 1, and larger M M makes the bound weaker. When increasing N, 
the lowerbound is naively expected to be scaled by \^N. Thus, the bound on M 3 from the 
Higgs boson mass is comparable to or rather stronger than that of the vacuum stability. 

Now we discuss the correlation between B(B — > X s -y) and for various Mm- It is 
shown in Fig. 3. The left (right) figure is for N = 1 (2). In these figures, the colored 
shaded regions are the experimental bounds at 1, 2, and 3 a levels, and the dark shaded 
region is the SM model prediction with a 1 a error. The red circle dots correspond to 
discrete M 3 's with 100 GeV intervals. The allowed regions are plotted with solid lines 
while the regions experimentally excluded by the light Higgs mass bound are also plotted 
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with dotted lines. 

We can see that the SM predictions for these processes are rather low compared to the 
experimental data. The sizes of the experimental errors are almost comparable to those for 
theoretical ones. Now we stand on a stage in which we can pin down a parameter region 
in models beyond the SM by combining experimental results of low-energy observables. 

As discussed in Section 3.1, the SUSY contribution to are always positive. The 
magnitude of the SUSY contribution is mainly determined by M 3 (A) and is not so sensitive 
to Mm and N. On the other hand, the SUSY contributions to B[B — > X s ^) depends on 
all the parameters. 

For smaller M M B(B — > X s ^y) is larger than the SM prediction, while B[B — > X s ^y) 
can be smaller than the SM prediction in relatively small M 3 region. We find that the 
chargino and the gluino contribution are not significant in small Mm regions, in which 
the non-SM contributions to C 7 are dominated by the charged Higgs contribution. Since 
the charged Higgs amplitude always contributes to C 7 constructively, B(B — > X s ^) in the 
MGM becomes larger in small Mm region. For large Mm region, the chargino contribution 
becomes significant and dominates the SUSY contribution to C-j. As a result, B(B — > X s ^) 
can be smaller than the SM prediction. 

We find that the most favorable region in the MGM from BiB — > X s '~f) and a M is 
M M ^ 2 A and M 3 ~ 900 (1200) GeV for N = 1 (2). The gluino mass will be covered by 
the direct SUSY search in the LHC experiment. The light Higgs boson mass is close to 
the current experimental lowerbound as in Fig. 2. 

In Fig. 4, we show the correlation between B(B — > X s ^) and B(B — > X s /i + /i"). The 
SUSY contributions to B{B — > X s fi + fi~) is negative and rather small. If the constraint 
from the light Higgs boson mass is taking into account, the deviation is at most few %. 
Considering both theoretical and experimental uncertainties, it would be very difficult to 
observe the deviation of B[B — > X s fi + fi~) in near future. 

In Fig. 5, we show the correlation between B(B — > X s '-f) and B(B S — > fi + fi~). B(B S — > 
becomes large for small Mm regions, since tan ft becomes very large (~ 40—50), and 
the tan 6 /3 enhancement on B(B S — > fi + fi~) is effective. When we consider the favorable 
region from B[B — > X s j) and a^, the branching ratio is estimated as B(B S — > /i + /i~) — 
(6 — 7) x 10~ 9 . This is still far below the current experimental bound, but with in a reach 
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of search at LHC. 

Finally, in Fig. 6, we show the correlation between B(B — > X s ^y) and B{B~ — > t~V). 
The deviation of B{B~ — > r~v) is always negative in the MGM and become significant 
for small M M regions. This is because tan/5 becomes very large in small M M and the 
charged Higgs contribution interferes destructively with the SM contribution. At the 
favorable region from B(B -> X s ^) and a„, B{B~ -^ri/)~lx 10~ 4 . Notice that this 
branching ratio is around the 1 a lower bound of current experimental data. A precious 
measurement of B{B~ — > r~ 17) at -B factories may be helpful for search for signature of 
the MGM. 

5 Summary and Discussion 

We have revisited phenomenology in the minimal gauge mediated model, in which tan/3 
is naturally large. We have considered the anomalous magnetic moment of the muon, 
the branching ratios of B — > X s j, B — > X s l + l~, B s — > and B~ — > t~V. When 

tan/? is very large, the SUSY contributions to those observables can be enhanced, and 
the deviations from the SM predictions are strongly correlated with each other. We 
have updated the results in [8] taking account of the recent theoretical and experimental 
developments of B — > X s ^ and the anomalous magnetic moment of the muon. 

We have shown that the experimental bound on the light Higgs boson gives a strong 
bound on the MGM and requires relatively heavy SUSY spectrum. We find that the 
lower bound of the gluino mass is M 3 > 800 ~ 900 (800 ~ 1100) GeV for JV = 1 (2). This 
constraint makes the signatures in the MGM less significant. 

When the messenger scale is larger than A by orders of magnitude, B[B — > X s ^y) is 
smaller than the SM prediction, since the chargino loop contribution interferes destruc- 
tively with the SM one. On the other hand, when the messenger scale is same order of A, 
B[B — > X s ^j) becomes larger than the SM prediction since the charged Higgs contribution 
becomes important. 

The most recent theoretical calculation of B(B — > X s ^y) in the SM [9, 10] is 1.4 a lower 
than the experimental world average [15]. When combining it with a 3.4 a deviation in a M , 
a parameter region in Mm ~ A and M3 ~ 1 TeV is favored. In the region, B{B~ — > t~ 77) 
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Figure 1: Two parameters, /x and tan/3, as functions of M 3 for N = 1 (left) and 2 (right). 
Each line in figures of [i parameter is for Mm = 2, 10, 10 2 , 10 3 , and 10 4 A from the bottom 
line. 
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M 3 [GeV] 
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Figure 2: Light neutral Higgs boson mass as a function of M 3 for N — 1 (left) and 2 
(right). Each line is for M M = 2, 10, 10 2 , 10 3 , and 10 4 A from the bottom line. 




Figure 3: Correlation between B{B — > X s 7) and a M in the MGM for N = 1 and 2. The 
colored shaded regions represent experimentally allowed ones at 1,2, and 3 a, respectively. 
The grey shaded region is theoretical prediction with 1 a uncertainty. The solid (dotted) 
lines are experimentally allowed (excluded) regions. The circle dots represent the discrete 
M 3 's with 100 GeV intervals. 
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Figure 6: Correlation between B(B — > X s , -f) and B(B — > r z/) in the MGM for N — 1 
and 2. 




Figure 7: Lightest stau and lightest neutralino masses as functions of M 3 for N = 1 (left) 
and 2 (right). 



16 



is predicted smaller than the SM prediction while B(B S — > is larger than the SM 

prediction. Precise measurements of B(B~ — > t~V) at B factories and B(B S — > /i + /U~) at 
TEVATRON and LHC may play important roles to search for signature of the MGM in 
addition to SUSY direct search at LHC. 

The nature of the next-lightest SUSY particle is important in the gauge-mediated 
models from viewpoints of collider physics and cosmology. In Fig. 7 we show the lightest 
stau and lightest neutralino masses as functions of M 3 for N — 1 and 2. The stau mass 
is lighter than or quite degenerate with the neutralino mass in the regions favored from 
B(B — > Xs'-f) and a^. This will have a important implication to SUSY particles searches 
and also the studies at LHC and LC. 

The parameter region M M ~ A is well-motivated from the cosmological gravitino 
problem. When the S field is directly coupled with the dynamical SUSY breaking sector 
and M M ~ 100 TeV, gravitino may be lighter than ~ 10 eV. Such an ultra-light gravitino 
is known harmless in cosmology [32]. 
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